THE FUJITA EXPONENT FOR SEMILINEAR HEAT 
EQUATIONS WITH QUADRATIC ALLY DECAYING 
POTENTIAL OR IN AN EXTERIOR DOMAIN 
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Abstract. Consider the equation 

u t = Ait - Vu + au p in R n x (0, T); 

(0.1) 

u(x, 0) = 0(a;) > 0, in R n , 
where p > 1, n > 2, T 6 (0, ool, Vfas) ~ i-*™ as lac I — > oo, for some w 7^ 0, 

|x| 

and a(a;) is on the order \x\ m as \x\ — > cxd, for some m G (—00,00). 
A solution to the above equation is called global if T = 00. Under 
some additional technical conditions, we calculate a critical exponent 
p* such that global solutions exist for p > p* , while for 1 < p < p* , 
all solutions blow up in finite time. We also show that when V = 0, 
the blow-up/global solution dichotomy for (|0.ip coincides with that for 
the corresponding problem in an exterior domain with the Dirichlet 
boundary condition, including the case in which p is equal to the critical 
exponent. 



1. Introduction and Statement of Results 
Consider the semilinear heat equation 

u t = Au - Vu + vP in R n x (0, T)\ 

(1.1) 

u(x,0) = 4>{x) > in R n , 

where p > 1, n > 1 and T G (0, cxd]. In this paper, when we speak of a 
solution to the above equation, or to any of the other equations appearing 
later on, we mean a classical solution u satisfying \ \u(-, t)\ |oo < oo, for 
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< t < T. This allows us to employ comparison principles. A solution to 
(jl.ip is called global if T = oo. In the case that V = 0, p* = 1 + ^ is the 
critical exponent, the so-called Fujita exponent, and one has the following 
dichotomy: if p > p* , then for sufficiently small initial data cj>, the solution 
to (jl.ip is global, whereas if 1 < p < p*, then (jl.ip has no global solution — 
every solution blows up in finite time. This result goes back to Fujita [3] in 
the case p ^ p* . Various proofs of blow-up in the borderline case p = p* can 
be found in PQ, [8], [H]. 

More recently, Zhang p3] considered (jl.ip with n > 3 for potentials V 
behaving like y^j^p, for 6 > and He proved the following result. 

Theorem (Zhang). Let n > 3. 

i. If < V(x) < /or some 6 > 2 and w > 0, iaen p* = 1 + | and 
consequently the potential does not affect the critical exponent; 

ii. If V{x) > y+j^p' for some b G (0, 2) and uj > 0, iaen p* = 1 and £aere 
exisi global solutions for all p > 1; 

Hi. If yqjfcfs < < 0, /or some 6 > 2 and a; < wrf/i |cj| sufficiently 

small, then p* = 1 + - and consequently the potential does not affect the 
critical exponent; 

iv. If V{x) < j-qq^p , /or some 6 G (0, 2) and a; < 0, then p* = oo and there 
are no global solutions for any p > 1. 

Note that wherever the statement of the result is that there exist global 
solutions, Zhang either does not allow for negative V or else requires that 
| V| be sufficiently small. The reason for this will become clear from Theorem 
[2 below. 

Zhang noted that it seemed difficult to specify the exact value of the 
critical exponent in the case of quadratic decay; that is in the case that 
V(x) ~ as \x\ — * oo. He also noted that it is unclear whether or not p* 
is finite in the case that V(x) ~ -pr?, with uj < 0. 
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Very recently, Ishige [61 treated (11. ip for n > 3 in the case V(x) ~ with 

\x\ 

lo > 0. Let a = a(u>,n) denote the larger root of the equation a(a + n — 2) = 
lu; that is 



(1.2) a(w,n) = ^ . 

Since we are assuming here that uj > 0, one has a(cj,n) > 0. Define 

(1.3) p» = l + — ^ T- 

n + a(w, n) 

Theorem (Ishige). Let n > 3 and assume that V > 0. Lei w > 0. 

i. IfV(x) > j^p /or /arge i/ien /or p > p*(u;) i/iere ezisi global solutions 
to jTTJ; 

ii. IfV(x) < ™ /or /an/e |x|, i/ien /or 1 < p < p*(tc>) every solution to 
(jl.ip blows up m finite time. 

Note that Ishige assumes from the outset that V > 0. The delicacy 
between having global solutions and allowing V to take negative values will 
be explained by Theorem [2] below. 

Ishige's proof involved comparison with a solution to the radially symmet- 
ric linear equation v% = Av — V(\x\)v, where V(r) ~ ^ as r —* oo. The large 
time behavior of this linear equation, which is needed for the comparison, 
was recently obtained by Ishige and Kawakami [TJ. 

In this paper, our main focus is the study of the remaining case, V(x) ~ 
r%j, with u < 0. In fact we treat the following more general problem: 

\x\ 

u t = Au-Vu + au p in R n x (0, T); 

(1.4) 

u(x,0) = 4>(x) > 0, in R n , 

where p > 1, n > 2, T € (0, oo], is bounded and continuous, < 
a 6 C a (R n ) and V 6 C a (R n - {0}) a G (0,1]. We also require that 
liminf 2 _i.o V(x) > — oo so that V is locally bounded from below. Our meth- 
ods, which are completely different from the method employed by Ishige, also 
allow one to obtain weaker versions of Ishige's results for the case to > 0, but 
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in the more general context of equation (jl.4[) with n > 2. The method of 
proof also leads naturally to a study of the critical exponent in an exterior 
domain with the Dirichlet boundary condition in the case V = 0. 

In the case that V = and that a satisfies 
(1.5) 

ci|x| m < a(x) < C2\x\ m , for sufficiently large \x\ and some m E (— oo,oo), c\,C2 > 0, 

the critical exponent p* for (|1.4[) was calculated in [12]; it is given by 

(2 + m)+ 



(1.6) p* = 1 + 



In (|1.2|) we defined a(w, n) for cj > 0. We now extend the definition of a(ui, n) 
in (TOD tow> -i(n - 2) 2 . Note that a(uj,n) < for -\(n - 2) 2 < w < 0. 
Now define 

n 7 , * / \ -i . (2 + m) + 

(1.7) p (to, m) = 1 + 



n + a(w, n) 
We will prove the following theorem. 

Theorem 1. Let n > 3 and Ze£ —\{n — 2) 2 < u; < 0. Consider (jl.4p wii/i 
a(x) satisfying fjl .5|) . Assume that V G C a {R n — {0}) and i/iai 
liminfa^o l^(x) > — oo. Letp*(uj,m) be as in (jl.7p . 

z. IfV(x) > |^p, iaen i/iere exzsi global solutions to (jl.4p /orp > p*(uj,m); 
ii. IfV(x) < T^p, /or sufficiently large \x\, then there are no global solutions 
to (II. 4p /or 1 < p < (cj, m) . 

Remark. Note that in the case of the existence of global solutions, we allow 
V to be negative up to a precise globally specified size. The reason for this 
will become clear in Theorem [2l 

We now consider what happens when u> < — l(n — 2) 2 , n > 2. We will 
show that p* = oo under a certain general condition on the operator —A+V, 
and that this condition holds if V(x) < j^p, for |x| > e, with sufficiently 
small e > 0. 
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Let D ^ R n be a domain. Then —A + V on D with the Dirichlet bound- 
ary condition on dD can be realized as a self-adjoint operator on L 2 (D). 
Denoting its spectrum by <r(— A + V;D), let 

X . D (-A + V) = mta{-A + V;D). 

Theorem 2. If there exists a domain D S R n for which ini x£ D a(x) > 
and \o-d(— A + V) < 0, then there are no global solutions to (jl,4p for any 
p > 1; that is, p* = oo. 

We can use Theorem [2] to prove the following corollary. 

Corollary 1. Consider (|1.4p with a > on R n , n > 2. Let uj < -\{n-2) 2 . 
There exists an e > such that ifV{x) < -0™, for \x\ > e, then there are no 
global solutions to (jl.4p for any p > 1; that is, p* = oo. 

Remark 1. Note that there is a discontinuity in the critical exponent at 

uj = —\{n — 2) 2 . By Theorem [Q if V(x) = — ^n^p , for sufficiently large 
and K(ic) > — > f° r an then the critical exponent is equal to 

j(n — 2) 2 ,m) = 1 + ■ However, if V(x) = -Sv, for some u> < 

— \{n — 2) 2 and \x\ > e, for sufficiently small e > 0, then the critical exponent 

is oo. 

Remark 2. Theorem [2] makes it clear why in the theorems of Zhang and of 
Ishige and in Theorem [TJ one needed to be careful with regard to stating 
the existence of global solutions and allowing V to take negative values. For 
example, part (iii) of the theorem of Zhang states that if ^nbs — V(x) < 
for some b > 2 and uj < 0, with \uj\ sufficiently small, then the critical 
exponent for (11.11) is 1 + \- The requirement that \uj\ be sufficiently small 
is mandatory in light of Theorem [2j Indeed, for any D Q R n , if uj < and 
\uj\ is sufficiently large, then Ao ; d(— A + j^nby ) < and thus, by Theorem 
[2J one has p* = oo. 

The method of proof in Theorem Q] also yields the following result for the 
case uj > 0. 
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Theorem 3. Let n > 2 and to > 0. Consider (II. 4|) with a{x) satisfying 
(fT3j) . Assume that V E C a (R n - {0}). Letp*(to,m) be as in (fLTj) . 
i IfV(x) > |^p, i/ten f/iere exisi global solutions to (jl.4p /orp > p*(u,m). 
ii. IfV(x) < T^p, /or sufficiently large \x\, then there are no global solutions 
to (PI) forl<p< p*(u),m). 

Remark. Note that part (i) requires that V approach oo as |x| — > 0. In fact, 
as Ishige has proven in the case that a = 1 and n > 3, the result should 
hold as long as V(a;) > ™, for sufficiently large \x\, and V > for all 
x. However, our method of proof does not seem to be extendable to this 
situation. 

As will be seen below, the method of proof we employ for the blow-up 
case in Theorems Q] and [3] will lead naturally to a consideration of the critical 
exponent for the semilinear heat equation in an exterior domain with the 
Dirichlet boundary condition and with V = 0. Let B r = {x £ R n : |x| < r}. 
Consider the following problem: 

u t = Au + au p in {R n - B ro ) x (0, T); 

(1.8) u(x,t) = 0, for \x\ = r , t > 0; 

u(x, 0) = 4>{x) > in R n - B ro , 

where 
(1.9) 

ci|x| m < a(x) < C2\x\ m , for sufficiently large \x\ and some m £ (— oo,oo), ci,C2 > 0. 

We prove that restricting to an exterior domain does not affect the blow- 
up/global solution dichotomy. 

Theorem 4. Let n > 2. Consider (jl.8p with a(x) satisfying f|l .9[) . Lei 

(2 + »)+ 

71 

as m (|1.6p . 

i- If ^ < P <P* , then there exist global solutions to (II. 8p : 
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it. If p > p* , then there are no global solutions to (jl .8() . 

Remark. In the case a = 1 , n > 3 and p ^ p* , the result in Theorem [4] 
was proven in [2]. For some other works that treat the critical exponent in 
exterior domains, see [9] and [15j . Most of the results in these papers do not 
cover the case in which p is equal to the critical exponent. 

We end this section with an outline of the methods used to prove Theo- 
rems Q] and [3l concentrating on the case of nonexistence of global solutions, 
which is where our method is novel, and leads to a consideration of the criti- 
cal exponent in the case of a semilinear heat equation in an exterior domain 
with the Dirichlet boundary condition. By standard comparison techniques, 
it suffices to treat the radially symmetric case. Thus, instead of considering 
solutions u(x,t) of ()1.4[) with a satisfying (jl.5p . we may consider solutions 
u(r, t) of the equation 

n — 1 

ut = u rr H u r — V{r)u + a{r)u p in (0, oo) x (0, T)\ 

(1.10) 

n(r,0) = <t>{r) > in [0,oo), 
where p > 1, T G (0, oo], <j) is bounded and continuous, V £ C a ((0, oo)) and 
liminf r _»o V(r) > — oo, < a £ C a ([0, oo)), a £ (0, 1], with a satisfying 
(1.11) 

c\r m < a{r) < C2r m , for sufficiently large r and some m £ (— oo, oo), ci, oi > 0. 

For the existence of global solutions when p > p*(uj, m) in part (i) of The- 
orems Q] and El we construct a global super-solution to (|1.1U|) . Note that in 
general it is much more difficult to use the method of super /sub-solutions to 
prove blow-up, since the construction of an appropriate sub-solution would 
probably require a reasonable knowledge of the blow-up profile. 

We now turn to the nonexistence of global solutions when 1 < p < 
p*(iv,m) in part (ii) of Theorems Q] and [31 We may assume without loss 
of generality that the initial data 4> m (jl-lOP satisfy 0(r) > for all r > 0. 
Indeed, if this is not the case, then for 5 > sufficiently small, we can con- 
sider u(r, t) = u(r,t + 5), which also satisfies (|1.10p and is strictly positive at 
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t = 0. We apply a transformation as follows. Let u be a solution to (jl.lOp 
and define v(r,t) = r~ a u(r,t). Let ip(r) = r _Q 0(r). Then one calculates 
that 
(1.12) 

v t = v rr + ^±^v r + ( a(a+ r r 2) - V(r)) - + r^a^ 

in (0,oo) x (0,T); 

t>(r, 0) = ip(r) > in (0, 00). 

There will be global solutions of v if and only if there are global solutions of 
u; thus it suffices to study (j!.12j) . In part (ii) of Theorems CD and El we are 
assuming that V(r) < for sufficiently large r, say for r > tq, where u > 
— \{n — 2) 2 . If one now chooses a = a(u>,n) as in (|1.2p . then the coefficient 
of v in (|1.12p is nonnegative for r > tq. By the comparison principle, the 
solution to that equation dominates the solution to the equation 

N — 1 



w t = w rr + 



-w r + a(r)w p in (tq, 00) x (0, T); 



(1.13) 



where 



w(r, 0) = ip(r) > in [7*0, 00); 
w(r ,t) = 0, i > 0, 



(1.14) N = n + 2a(uj,n) 

and a(r) = r a<yP ~ 1 ^ a{r) . In terms of a, the assumption (11 . 1 If) on a is 

(1.15) c\r M < a(r) < c<ir M , for sufficiently large r, ci,C2 > 0, 
where 



(1.16) 



M = a(u, n)(p — 1) + m, m 6 (—00, 00). 



(The reason we insisted on ^(r) > for all r > 0, and thus also ip(r) > for 
all r > 0, is that otherwise we could have ended up with tp = in (II. 13ft .) 
Thus, it suffices to show that there are no global solutions to (|1.13p - (|1.16p . 
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Now (fTT^ - ffTTTBl ) is the radial version of ([I3J-CL5]) in the case V = 0, 
except that we have placed the Dirichlet boundary condition at r = ro 
instead of considering the problem for all r > 0, and except that m in (jl.5jl 
is replaced by M and the dimension n is replaced by the "dimension" N. 
(Note from the definition of a(co, n) that one always has N > 2.) The critical 
exponent p* for (jl.4p with V = and with a satisfying (jl.5p was given in 
(jl.6p . Substituting iV and M for n and m in (|1.6p . it is not unreasonable to 
suspect that no global solutions will exist if 

(1.17) 1 < P < 1 + (2 + M)+ = 1 + (2 + ^ n)(P " 1} + m)+ 



JV ra + 2a(w,n) 

We now solve (j!.17j) for p. Consider first the case oj < 0, in which case 
a(uj,n) < 0. Since we are assuming that p > 1, (|l,17p will never hold if 
2 + a(u>,n)(p — 1) + m < 0; that is, if 

2 + m 

(1.18) P>l--7 r- 

a{u), n) 

On the other hand, if 2 + a(u>, n)(p — 1) + m > 0, then solving (j!.17jl for p 
gives 

2 + m 

1.19 Kp<l+ -, r. 

n + a(ui, n) 

One can check that for m > — 2, the right hand side of (|1.19p is strictly less 
than the right hand side of (|1.18|) . From this fact along with f|l. 18j) and 
(|1.19p . we conclude that (|1.17p holds if and only if 1 < p < p*(u, m), where 
p*(uj,m) is as in (|1.7|) . 

Now consider the case to > 0, in which case a(ui, n) > 0. 
If 2 + a(uj,n)(p — 1) + m > 0, then solving (|1.17|) as we did above gives 
(|1.19|) . On the other hand, if 2 + a(uj, n){p — 1) + m < (which implies 
that m < —2), then (|1.17p does not hold. Putting these facts together leads 
again to (|1.17p holding if and only if 1 < p < p*(co, m). 

To turn the above argument into a rigorous proof, we need to show that 
indeed no global solutions exist for (|1.13p - (|1.16p when 1 < p < 1 + ^ 2+ ^ ■ 
That is we need to show that the proof in |12J, which treated the operator A 
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in R n (whose radial part is 4p + can accommodate two changes: (1) 

operators of the form + ^7-^ with fractional N and (2) the Dirichlet 
boundary condition at r = tq, which serves to make solutions smaller. The 
proof in [12] made rather heavy use of the explicit form of the heat kernel 
p(t,x,y) = (47rt) _ 2" exp(— ) for the corresponding linear operator A — 
01 in R n . In the present case, the corresponding linear operator is + 
^~F^"§F ~ M wr th the Dirichlet boundary condition at r = tq. It turns out 
that if N > 2 (equivalently, u> > — \{n — 2) 2 ), then the heat kernel for 
this operator is comparable in an appropriate sense to the heat kernel for 
J^j + ^jr--§^ — -§1 on the entire space r > 0; thus, we will be able to use 
this latter heat kernel, which we can exhibit explicitly. However, this latter 
heat kernel is a much less convenient object than the Gaussian heat kernel. 
In fact, this obstacle prevented us from using the method of proof in [12] to 
prove the existence of global solutions above the critical exponent; hence the 
use of super-solutions. However, we were able to use this heat kernel and 
amend the nonexistence proof in [12] at or below the critical value. When 
N = 2 (equivalently, u> = —\{n — 2) 2 ), the heat kernel with the Dirichlet 
boundary condition is not comparable to the heat kernel on the whole space, 
however an appropriate lower bound is known and sufficient for our needs. 

In section 2 we prove the existence of global solutions in part (i) of The- 
orems [H and EJ In section 3 we prove the nonexistence of global solutions in 
part (ii) of Theorems [T] and [31 In section 4 we prove Theorem [U In section 
5 we prove Theorem [2] and Corollary [TJ 

2. Proofs of Part (i) of Theorems [U and [3] 

We assume that p > p*(uj,m), where p*(uj,m) is as in (|1.7|) . As noted 
in the first section of the paper, instead of studying (jl.4p with a satisfying 
(|1.5p . it suffices to study the radial problem (jl.lOp with a satisfying (jl.lip . 
By the standard theory, it suffices to exhibit a global super-solution. We 
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look for such a super-solution in the form 



n 1 

r cr ^ 



v(r, t) = 5- -t— exp( -), 

v ' ; (t + 1)1 y t + 1 

for some 5, c > and some a,j € (— oo, oo). We have 

,a 2cr . 

(2.1) v r = (- - — )«; 



, .a 2 4c 2 r 2 4ca a 2c . 

(2 - 2) Vrr = ^ + J¥TW~T + -i--*-—i )v] 



7 cr 2 

(2 - 3 > "' = <-*tt + (Tm)>' 

The condition on V in part (i) of Theorems Q] and [3] is that V(r) > with 
_^n-2)_ < w < g j n Theorem [T] and u> > in Theorem Using this along 
with (|2"3]h d22} and (f%3jl . we have 

^ 

ti _1 (r; r , r -| v r — V(r)v — vt + a(r)u p ) < 

r 

. , 9 . r 2 a 2 + (n — 2)a — uj 7 — 4ca — 2cn 

(2.4) ^ - + — — + rn 

+ ^Vr), ^ n exp(-^- 1 ' r2 ). 

V ; (t + l)7(P~l) PV t + 1 

In order to make the first term on the right hand side of (|2.4p vanish, we 
choose c = j, and in order to make the second term on the right hand side 
of (|2.4p vanish, we choose a = a(u>, n) as in (|1 .2|) . 

If m < 0, the assumption on a in (jl.lip guarantees that for some C > 0, 
a(r) < Cr m , for all r > 0. If m > 0, the assumption on a in (jl.lip guarantees 
that for some C > 0, a(r) < C(r V 1) TO , for all r > 0. This forces us to break 
up the next part of the proof into two cases. We will continue the proof 
under the assumption that m < 0. After the completion of this case, it will 
be easy to point out how to handle the case m > 0. 

Since a{r) < Cr m , the final term on the right hand side of (|2.4p (with c = 
I and a = a(u>, n)) is bounded from above by Cd^ 1 ^^^^^- exp(— ] 
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2 

Letting z = this upper bound can be written as 

(t _|_ l)(7-|«(w,n))(p-l)-|m 

which is itself bounded from above by ; — fi^ ? , w — ; — i — , where 

Ci = sup 2>0 -2 i a (^. n )(p- 1 )+| m exp(— I (p — l)z). In light of the above analy- 
sis, it follows from (|2.4p that 



„a(u>,n) r 2 
(2-5) v(r,t)=5———eM 



satisfies 



(2.6) 



If 



(t+l)i " v 4(t + l) 



v 1 (w rr H fr — V(r)f — v t + a(r)u p ) < 

7-a(w,n)-|n CxCSP- 1 



t+1 U+l\('Y-2 a (. u ',n))(p-i)-%'i 



(2.7) 7 - a(w,n) - < 
and 

(2.8) ( 7 -ia(u;,n))(p-l)-im>l, 

then after choosing 5 > sufficiently small, the right hand side of (|2.6p will 
be non-positive. The two inequalities (|2.7p and (|2.8p together are equivalent 
to 

1 1 -j_ I m ]_ 
-a(u,n) H ^— < 7 < a(u>,n) + -n, 

2 p — 1 2 

and this latter pair of inequalities can be solved for 7 if and only if ot(u, n) + 
< 2a(w,n) + n, or equivalently, if and only if p > 1 + n+ 2 ^™ w ) ■ Since 
we have assumed from the outset that p > 1, we conclude that if p > 
1 + n ) = then it is possible to choose 7 so that (|2.7p and 

(EBD hold. 



In the case m > 0, we have a(r) < Cr m , if r > 1, and a(r) < Cr°, if 
< r < 1. Thus, in order for the above analysis to go through in this case, 
we need to have (|2.8p hold as it is written and also with m replaced by 0. 
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However, since m > 0, if (12. 8p holds as it is written, then it holds a fortiori 
with m replaced by 0. 

In the case u> > 0, the function v given by (12. 5p with 5 > sufficiently 
small and 7 chosen to satisfy (|2.7[) and (|2.8p serves as an appropriate global 
super-solution. 

In the case u < 0, there is one technical problem; namely, that a(u, n) < 
and thus v is not finite at r = 0. This artificial singularity arises from the 
use of polar coordinates. Unfortunately, if one replaces r by r + c for some 
c > 0, then u will no longer be a super-solution. Thus, we argue as follows. 
Consider u and p > p*{oo, n) as fixed. Our work so far allows us to conclude 
that for sufficiently small initial data <j>, the solution u(x, t) of fll-4|) satisfies 
u(x,t) < v(\x\,t) up until some possibly finite blow-up time. Choose e > 
sufficiently small so that p > p*{uj — e, n). The function v in (|2.5p was shown 
to be a super-solution for (jl.4p under the assumption that the potential 

V satisfies V(a;) > i^p. Recall that in (|1.4|) we are also assuming that 

V is locally bounded from below. Therefore, there exists an ro > such 
that V(x) > for \x\ < tq. One can check that it is then possible to 
choose an xo / such that V(x) > , zr^n ■ Now consider the radial version 
(If .1 0p of p.4p but with the origin shifted to the point Xq. Call the new 
radial variable p = \x — xq\. Since we have V(p) > the construction 
above shows that there exists a function v(p,t) = ^ p ^ t+1 ^ ex P(~ 4(t+i) ) 
such that for sufficiently small initial data cj>, the solution u(x, t) of (jl.4p 
satisfies u(x,t) < v(\x — xo\,t) up until some possibly finite blow-up time. 
We conclude that for sufficiently small initial data <ft, the solution u(x, t) 
of (|1.4j) satisfies u(x,t) < v(\x — xo\,t) A v(|x|,t) up until its blow-up time. 
But the right hand side is finite for all x and t. Thus u is in fact a global 
solution. 
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3. Proofs of Part (ii) of Theorems [Q and [3] 

As was shown at the end of section 1, in order to prove that when 1 < 
P < p*(uj,m) there are no global solutions to (jl.4|) with a satisfying (jl.5p . 
it suffices to show that there are no global solutions for (|1 . 13j) - ()l . 16j) when 
p satisfies (|1.17j) . We will always assume that M > — 2 since otherwise 
there is nothing to prove. We wish to employ the method of proof used in 
|12j . This method requires a fairly explicit knowledge of the heat kernel for 
the corresponding linear equation. In the present case, the linear equation 
is Wt = W rr + ^-^-W r with (r, t) € (ro,oo) x (0, oo), for some possibly 
fractional N with N > 2, and with the Dirichlet boundary condition at 
r = ro. Denote the heat kernel for this equation by q~(N,r )(t-, r i P)- 

Denote by qm\(t, r, p) the heat kernel for the equation Wt = W rr + ^^-W r 
with (r,t) S (0, oo) x (0, oo). The kernel q(N){t, r i p) is the transition prob- 
ability density for the Bessel process of order N, and is given by [5] 

r 2 + p 2 p N ~^ rp 

(3.1) q( N ){t,r,p) =exp( — - ) fy ( ) 

4t 2t{rp)-~ 1 2 2t 

where l v is the modified Bessel function of order v , given by 

oo 

(3.2) i v {x) = {\yy j 



n=0 



n\T{u + n + 1) 



By the maximum principle, q^^(t,r, p) < q(N)(t, r, p). What we need, 
however, is an appropriate inequality in the reverse direction. 

If N > 2 (equivalently, uj > —\{n — 2) 2 ), then the Bessel process corre- 
sponding to the operator + ^jr--^ is transient [TT]. Furthermore, as will 
be explained momentarily, the uniform parabolic Harnack inequality holds 
for the heat equation Wt = W rr + "T 1 W r on r > 0. Thus, it follows from 
[3] that there exist constants Kq , c > such that 
(3.3) 

<?(7V,ro)(^ r >P) > cQ(N)( K ot,r,p), for r > r + l, p > r + l, t > and iV > 2. 
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(The uniform parabolic Harnack inequality concerns nonnegative solutions 
W of W t = W rr + W r on r > on a time interval [r, t + T) . See gl 
Definition 2.2] for the precise definition. Any such solution can be repre- 
sented as W(r,r + t) = J* °° qrm(t,r, p)W(p,r)dp, < t < T. Using the 
explicit formula for in (|3.1[) . one can verify the uniform parabolic Har- 
nack inequality. Indeed, in the case that N is an integer, the above heat 
equation is just the radial form of the standard heat equation on R N , and it 
is well-known that the uniform Harnack inequality holds in this case |10|.) 

The following key a priori lower bound on solutions to (|1.13|) - (jl.l6p in the 
case that 1 < p < 1 + - — N will be used to prove the theorem. Then we 
will come back to prove the lemma. 

Lemma 1. Let w be a solution to f|1.13j) - f|1.16() on a time interval < t < T, 

fM 
N 



with \<p<\ + (2+ ff } and N > 2. Then for some K, C > 0, 



n Kr 2 
(3.4) w(r, t) > Cr~ log(l + t) exp( — ), for 2 < t < T, r > r + 1. 

Remark. The proof of Lemma [1] makes use of (|3.3p . If N = 2, a weaker 
lower bound holds for q~(N,r ) i n terms of gYjv)- This weaker bound is enough 
to prove (|3.4p when N = 2 with the restriction that r > t? . See Lemma [2] 
and (|4.2p . As the proof of Theorem Q] below shows, it is enough to have the 
estimate (|3.4p for r > ia. 

In light of the above remark, (|3.4j) holds for all iV > 2 and r > t%. We 
now use this to prove the theorem. 

Proof of Theorem [0 Assume that w(r, t) is a global solution to (|1.13j) - (|1.16j) . 
For n > rp + 1, define 

F n {t)= w{r,t)^ n) {r)r N - l dr, 

J n 

where > 0, normalized by f^ n ^ n \r)r N ~ l dr = 1, is the eigenfunction 

i d 2 



corresponding to the principal eigenvalue A n > for the operator + 



^T^~1F-) = ~ rl N ~§r rN l ~dlr 011 ( n ' w itk the Dirichlet boundary condition 
at the endpoints. For an appropriate value of n, we will show that F n blows 



16 ROSS PINSKY 

up in finite time, thereby contradicting the assumption that w is a global 
solution. 

From the outset, we assume that n is sufficiently large so that (|1.15j) 
holds for r > n. Simple scaling shows that A n is on the order 4j as n — > 
oo. In particular then, there exists a constant c > such that A n < 
Since cp^(n) = ^ n \2n) = 0, one has ((p^)'(n) > and {^)'{2n) < 0. 
Using the facts in this paragraph, integrating by parts and using Jensen's 
inequality, we have 



<-2n 

F' n (t)= / wt{r,t)<j>^{r)r^- x dr 



n 



(r,t) + — - 1 w r (r,t) + a(r)w p (r,t)^J <p {n) (r)r N ~ l dr 



n 

2n r 2n 



(3.5) 



(r N - 1 w r (r,t)) r (t ) {n) (r)dr+ / d(r)w p (r, t)^ n) {r^^dr 



>/ {r N - x $f\r)\w{r,t)dr + c x n M \ w v {r ,t)4 n \r)r N ~ l dr 



r-2n 

-X n F n (t) + Cl n M / w p (r, t)^ (r^^dr 



>-—F n (t) + Cl n M FP(t). 



n 2 



The function — -^x + c\n M x v is both positive and increasing for x > 

1 M+2 

(c~) p ln p 1 • Therefore, if there exists an n and a T n for which F n (T n ) > 

1 M+2 

(■^) p ~ 1 n p- 1 , then it follows from (|3.5p and the fact that p > 1 that F n (t) 
will blow up at some finite value of t. From Lemma Q] and the remark 
following it, we obtain w(r,n 2 ) > C\n~ N \ogn, for n < r < 2n and some 
Ci > 0. Thus, F n {n 2 ) > C x n- N \ogn. Since 1 < p < 1 + (2+ ^ /)+ , one can 

-N M+2 

choose n sufficiently large so that F n (n ) > C\n log ra > (^-) p_1 n f- 1 . 

□ 



Proof of LemmaUl The solution W to the corresponding linear problem 
Wt = W rr + ^pW r with the Dirichlet boundary condition at r = ro and 



FUJITA EXPONENT FOR SEMILINEAR HEAT EQUATION 17 

with initial data tp is given by 

(3.6) W(r,t)= q( N>ro )(t,r,p)il}(p)dp. 

J ro 

By comparison, the solution w to (|1.13j) - (|1.16p satisfies 

(3.7) w > W. 

On the other hand, the solution w to (|1.13p - (|1.16p satisfies the inequality 
(3.8) 

roo f*t /*oo 

w(r,t)> q(N,r )(t,r,p)^(p)dp+ ds dp q {N ^ o) (t-s,r, p)d(p)w p (p, s). 
Jvq J J ro 

(See [12] and [13], where it is also shown that under appropriate conditions, 
(|3.8p holds with an equality.) Without loss of generality, we assume that 
ro + 2 is contained in the support of ip appearing in (|3.6p . From (|3.ip - (|3,3p 
and (f376|) - (f378|) it then follows that 
(3.9) 

f't f'OO 

w(r,t) > ci / ds dp q {N) {Ki(t-s),r, p)d{p)q p {N) {K l s, p,r G +2), r > r +l, 

JO Jro+l 

for some K\,c\ > 0. 

In the case that N is an integer, which we denote by Nq, c/(at ) is just the 
standard iVo-dimensional Gaussian heat kernel in radial coordinates, and 
(|3.9p can be converted back to iVo-dimensional Euclidean coordinates. In 
[12], the right hand side of (|3.9p (converted to Euclidean coordinates and 
with some other inessential differences) was shown to satisfy the inequality 



(3.10) 



ds I dp q( No) (Ki(t - s),r,p)d(p)q p {N JK 1 s,p,r + 2) > 

ro+1 

.!*0. „j.M , Kr .2^ .„ _ 2+M 



Ct 1 - 2 p+- exp(-£p), if p< 1 + 



N 



Ct-^]og(l + t)exp(-Z£), ifp=l + 2±M, 



for i > 2, r > ro + 1, 

where K,C > 0. Recall that we are assuming that M > —2. Note that 
1 _ ^ p + f > _^Vo , if p < 1 + and 1 - + f = -fi, if p = 

1 + TvT- Thus ' from and <I3TH it follows immediately that dSHJ) holds 
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for N = N . (For ff37TO|) and (pT4"|) with N = N , see the statements and 
proofs of [12 j Lemma 2, Proposition 1 and Lemma 3]. The spatial integral 
in [12] is over all of R N °, which would correspond here to p > 0. But one 
could have worked just as well with |x| > ro + 1 in [12], so the restriction 
here to p > ro + 1 in the spatial integral causes no problem.) 

We now proceed to demonstrate that (|3.10p . and consequently also (|3.4p . 
continue to hold in the case that Nq is replaced by any non-integral N > 2. 
We write N = Nq — [3, where iVo > 3 is an integer and f3 £ (0, 1). Let 
K v {x) = I u (x), and note from the definition of I u in (|3.2p that K v (x) 
is decreasing in v. Thus, we have from (|3.ip . 
(3.11) 

u \ i r + P \ P t r P\^—\js 
1(N)(t,r,p) = exp( — - ) - (— ) 2 -fTiv 

4t 2t(rp)~- 1 2t 
r 2 +p 2 p^- 1 f rp rp 

> ^ ft ^ 

> — a-<l(N ){t,r,p). 

From (|3.1ip we have 

Q(N)( K l(t ~ s),r, p)q p {N) {K 1 s,p,r + 2) 

(3-12) > C 1 t ^s%\ No) (K 1 (t-s),r,p)<f m (K 1 s,p,r + 2), 

for 1 < s < — , 0<p< oo, 
- - 2 ~ 



t 

2 

ds 




for some Ci > 0. From (|4.8p we have 
(3.13) 

rfp Q(N)( K i(t ~ s),r, p)a(p)q p (N) (Kis, p,r + 2) > 



ro+l 



Cita / ds dp p P S 2 p q {No) {Ki{t- s),r,p)a(p)q P No) (K 1 s,p,r + 2). 

Jl Jro+l 

Note that the only difference between the terms appearing inside the 
double integral on the right hand side of (|3.13p and the terms appearing 
inside the double integral on the left hand side of (I3.10P is the addition of 



FUJITA EXPONENT FOR SEMILINEAR HEAT EQUATION 19 
— ft S- 

the factors p p and s 2 . Translating the setup and notation in the proof of 
(|3.1(J|) in [12J to the present situation, we note that the integration over p 
introduced a term of the form ((i — s)r(s, t))~z , where r(s,t) = s+pK s 2 ^ t _ s ^ , 
for some K2 > 0, and the exponent 4^- was a consequence of a being on 
the order p M . Since a(p) is replaced by p~^a{p) in (|3.13p . in the present 
situation we obtain a term of the form ((i — s)r(s, t))~~ 2; see |12^ (2.34)- 



(2.37)]. Thus, whereas in the penultimate step in the proof of (|3.10p in [12 
we obtained 



ds I dp q( No )(Kx(t- s),r,p)a(p)<fl N JK 1 s,p,r + 2)> 

Jro+l 



M 



C 2 exp( ) / s~— p (r(s,t))— + ^(t-s)^ds, 

for some K > (see [12, (2.37)]), we obtain here 
(3.14) 

f~2 f°° /3 

ds dp p- p s* p q( No) {Kx(t- s),r, p)d(p)q P No) (K lS , p,r + 2) 
Jl Jro+l 

>C 2 * 2 exp( ) / s 2 p ^2P[r(s,t)) 2 + 2 2 (t — s) 2 2ds. 

t Ji 

Making the change of variables u = § and recalling that Nq — (3 = N, we 
have 

if / s~^ Ip+ § p (r(s, t))^"" 1- ^ - 2 (i — s)^~2ds = 



(3.15) 



1 

1 1 iVT _ iV /" 2 JV,M_JV 
t^2 / U2+2 2P(tJ +pil2(l — U)) 2 2 (1 — u) 2 2 cfoj. 



If p < 1+ 2 "^- , then -f- + 4f — > — 1 and the integral on the right hand side 
of (ETT5]) is bounded in t. However if p = 1 + 2±M ; then f + 4f - f p = -1 
and that integral is on the order of logt. Using this fact along with ()3. 13[) - 
(I3.15p . we conclude that (13. 1QH holds with the integer Nq replaced by non- 
integral N. From this and f|3.9[) we then also obtain (|3.4p with the integer 
Nq replaced by non-integral N. This completes the proof of Lemma [TJ 

□ 
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4. Proof of Theorem 0] 

Note that (|1.13p - (|1.16|) with N equal to an integer is the radial version of 
(jl.8p - (|1.9j) (with N and M identified with n and to). Thus, in fact, Lemma 
[T] and the proof of Theorem [1] given in section 3 give a proof of Theorem [J] 
in the case n > 3. If we prove the equivalent of Lemma [1] for n = 2, then 
we will also have a proof of Theorem |4] for n = 2. In fact, as the proof of 
Theorem Q] showed, it suffices to have the estimate on w(r,t) in Lemma Q] 
for r > ia. Thus, it suffices to prove the following result. 



Lemma 2. Let w be a solution to (jl.8p «ra£/t n = 2 on a time interval 
< t < T, with 1< p < 1 + (2+ " )+ . T/ien /or some K, C > 0, 

i^lrrl 2 i 
(4.1) u;(x, t) > Cr l log(l + t) exp( L-L), /or \x\ > 1 2 and5 <t <T. 

Proof. We assume that m > —2 since otherwise there is nothing to prove. 
Let p(t,x,y) = (47rt)~ 1 exp(— ^ ) denote the heat kernel for the Lapla- 
cian on R 2 , and let p ro (t,x,y) denote the corresponding heat kernel for the 
Laplacian on R 2 — B ro with the Dirichlet boundary condition at \x\ = tq. It 
was shown in [3] that for appropriate constants cq, Kq > 0, one has 
(4.2) 

« ( txv )> cn log(l + N)log(l + ||/|) 

W ' °(log(l + ^) + log(l + N))(log(l + ^)+log(l + |y|) n °' 

for |x| > ro + 1, \y\ > ro + 1, i > 0. 



We now follow to a significant degree the proof of blow-up in [12]. Similar 
to |12[ Lemma 1], we have 
(4.3) 

w(x,t)> / Pr (t,x,y)(j)(y)dy+ / p ro (t-s,x,y)a(y)w p (y,s)dyds 



ro V) yj K r\yj u 'yi / / Fro v 

R 2 -B rQ JO JR 2 -B r , 



The first term on the right hand side of (|4.3p . which is the solution of the 
corresponding linear problem, constitutes a lower bound for w. Thus, using 
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(|4.2p . we have similar to p2J Lemma 2], 

1 / \x\ 2 -> log(l + \x\) 

w{x,t) > ct _1 exp(- ! ' ^ 



2if r (l og (l + Vt) + log(l + \x\)) (log(l + Vt)) 

for some c > 0. Note that for \x\ > £2, |y| > and £ > 1, the expression 
t, — - — f~. r^r-^vr^%j-^7r; , — n — rrv is bounded and bounded away from 

(log(l+v / t)+log(l+|x|))(log(l+ v / t)+log(l+|y|) J 

0. Thus, substituting the estimate (|4.4p into the second term on the right 
hand side of (|4.3I) . and using (|4.2p and ()1.9p . it follows that for some C > 0, 

c f ,™ , b-z| 2 « / |y|V 



w(x,t)>— / s p |y| m exp( — — ) exp(-— - )dyds, 

(4.5) 1 J\y\>ti ct 2K o s 



i 



for \x\ > £2 and large t. 

Performing some algebraic manipulations similar to those in [121 p. 166], one 
has for £, s > 1 and some c > 0, 

I 1 i 1 2 i 1 2 i 1 2 

y — a; y p. m y 

(4.6) exp( — — ) exp(--— ) > exp( -) exp( ). 

Ct AKqs ct cs 

Recalling that m > —2 and that n = 2, it is not hard to show, similar to 
Lemma 4], that for some k > 0, 

|2 

i i m ( \y 
j 1 2/| exp( 

!\y\>t- 

From (|4.5p - (|4.7p . we obtain for some k\ > 0, 



(4.7) / . |y| m exp(-^)dy > ks 1+ ^ , for s > ii 

J\v\>ti cs 



(4.8) w(x,t) > — exp(— — — ) / s 1+ ^"~ p ds, for |x| > t a and large £. 

By assumption, 1 < p < 1 + 2±p = 2 + f ; thus, 1 + f - p > -1. Conse- 
quently, for some &2 > and t > 5, we have 

(4.9) / 1 s 1+ ^~ p ds > k 2 logt. 

Jt2 

Now (|4TT|) follows from P~8j) and (JOJ. □ 
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5. Proofs of Theorem [2] and Corollary □ 

Proof of Theorem^ It is known that Xq-d(~ A + V) is non-increasing in D 
and that A 0; d(-A + V) = lim^^ A 0; D fc (-A + V), if L> fc | D p3J chapter 
4]. These properties of Ao ; d(— A + V) allow us to assume without loss of 
generality that the domain D in the statement of the theorem is bounded and 
has a smooth boundary. As such, Ao ; d( — A + V) < is in fact the principal 
eigenvalue for —A+V in D with the Dirichlet boundary condition. Let ipo > 
0, normalized by J D ^o(x)dx = 1, denote the corresponding eigenfunction. 

Assume now that u[r, t) is a global solution to (|1.4|) for some p > 1. 
Define 

(5.1) = / u(x,t)ip (x)dx. 

Jd 

We will show that F blows up at some finite time, thereby contradicting the 
assumption that u is a global solution. Note that ipo vanishes on dD and 
that VV'o • v < on dD, where v is the unit outward normal to D at dD. 
Also, by assumption inf^go a(x) > 5, for some 5 > 0. Integrating by parts, 
and using Jensen's inequality and the facts above, we have 

F'(t)= / Ut(x, t)ipo(x)dx = / (Au — Vv + au p )(x)^o(x)dx 

(5.2) Jd Jd 

> -A 0; d(-A + V)F(t) + 5F p (t) > 5F p (t). 

Although the initial data <j) of u may vanish identically on D, one certainly 
has F(t) > for t > 0. Thus, it follows from (|5.2p and the fact that p > 1 
that i* 1 blows up at some finite time. □ 

Proof of Corollary [3 It is well-known that A 0; /jn_{o}(— A + j^j) < if 
7 > ^r^- HI PP- 153-154]. Let B k = {x G R n : \x\ < k}. Recalling 
the facts noted in the first line of the proof of Theorem [21 it follows that 
A 0; B fc _B E (— A + jX;) < 0, for sufficiently large k and sufficiently small e > 0. 
Since a is continuous and positive by assumption, it follows that a is bounded 
away from on B k — B e . Thus, the corollary follows from Theorem [2j □ 
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